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Abstract 

Let G be a graph with adjacency matrix A(G). Consider the matrix 
IA(G) = (I | A(G)), where I is the identity matrix, and let M(IA(G)) 
be the binary matroid represented by IA(G). Then suitable parametrized 
versions of the Tutte polynomial of M(IA(G)) yield the interlace polyno- 
mials of G, introduced by Arratia, Bollobas and Sorkin [J. Combin. The- 
ory Ser. B 92 (2004) 199-233; Combinatorica 24 (2004) 567-584]. Interlace 
polynomials subsequently introduced by other authors may be obtained 
from parametrized Tutte polynomials of the binary matroid represented 
by (/ | A(G) | I + A(G)). 

Keywords, interlace polynomial, matroid, multimatroid, Tutte poly- 
nomial 
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1 Introduction 

Motivated by problems that arise in the study of DNA sequencing, Arratia, 
Bollobas and Sorkin introduced a one-variable graph polynomial, the vertex- 
nullity interlace polynomial, in [5] . In subsequent work [3J 3] they observed that 
this one- variable polynomial may be obtained from the Tutte-Martin polynomial 
of isotropic systems studied by Bouchet [8] , introduced an extended two- variable 
version of the interlace polynomial, and observed that the interlace polynomials 
are given by formulas that involve the nullities of matrices over the two-element 
field, GF(2). Inspired by these ideas, Aigner and van der Hoist [JJ, Courcelle [13] 
and the author [551 US] introduced several different variations on the interlace 
polynomial theme. 

All these references share the underlying presumption that although the 
theory of the interlace polynomials is connected to that of the Tutte polynomial 
in some ways, the two theories are largely separate in general. In this short note 
we point out that in fact, the interlace polynomials of graphs can be derived from 
parametrized Tutte polynomials of binary matroids associated with adjacency 
matrices. We presume the reader is familiar with the standard terminology of 
graph theory and matroid theory; see [TBI 22, 28, 29] for instance. 
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2 The identity-adjacency matroid 



We restrict our attention to looped simple graphs. That is, a graph G is given by 
specifying a finite set V(G) of vertices, declaring that certain vertices are looped 
and the others are not, and declaring that certain pairs of distinct vertices 
are neighbors and the other pairs are not. The adjacency matrix of G is the 
V(G) x V(G) matrix A(G) with entries in GF(2) given by: a diagonal entry is 
1 if and only if the corresponding vertex is looped, and an off-diagonal entry is 
1 if and only if the corresponding vertices are neighbors. 



If 7 is the \V(G)\x \V(G)\ identity matrix then IA(G) is the \V(G)\x2 \V{G)\ 



For convenience of notation and in order to indicate the relationship with our 
previous work on the interlace polynomials, we use the Greek letters <j> and \ 
to refer to the columns of the indicated submatrices of IA(G): the column of I 
corresponding to v is denoted ity, and the column of A(G) corresponding to v 
is denoted v x . 

Definition 1 The identity-adjacency matroid M(IA(G)) is the binary matroid 
represented by IA(G). 

That is, M(IA(G)) is a matroid on the ground set W(G) = {v^, v x | v € 
V(G)}, and if T C W(G) then the rank r G {T) of T in M(IA{G)) equals the 
dimension of the GF(2)-vector space spanned by the columns of IA(G) corre- 
sponding to elements of T. 

One way to define the Tutte polynomial of M(IA(G)) is a polynomial in the 
variables s and z, given by the subset expansion 



We do not give a general account of this famous invariant of graphs and matroids 
here; thorough introductions may be found in [5J Q21 [TH [TB] . 

Tutte polynomials of graphs and matroids are remarkable both for the amount 
of structural information they contain and for the range of applications in which 
they appear. Some applications (electrical circuits, knot theory, network relia- 
bility, and statistical mechanics, for instance) involve graphs or networks whose 
vertices or edges have special attributes of some kind - impedances and resis- 
tances in circuits, crossing types in knot diagrams, probabilities of failure and 
successful operation in reliability, bond strengths in statistical mechanics. A 
natural way to think of these attributes is to allow each element to carry two 
parameters, a and b say, with a contributing to the terms of the Tutte polyno- 
mial corresponding to subsets that include the given element, and b contributing 
to the terms of the Tutte polynomial corresponding to subsets that do not. Za- 
slavsky [30] calls the resulting polynomial 



matrix 



IA(G) = (I | A(G)) . 



t(M(IA(G))= J2 s rGWG)) - rG(T) z |T| - rG(T) . 



TCW(G) 



TCW(G) 




T G {W(G))-r G {T) z \T\-r G (T) 



(1) 
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the parametrized rank polynomial of M(IA(G)); we denote it t(M(IA(G)). 

We do not give a general account of the theory of parametrized Tutte poly- 
nomials here; the interested reader is referred to the literature, for instance 
[7J [TBI [531 (Ml ISO]- However it is worth taking a moment to observe that 
parametrized polynomials are very flexible, and the same information can be 
formulated in many ways. For instance if s and the parameter values b(w) are 
all invertible then formula (flj is equivalent to 



which expresses t(M(IA(G)) as the product of a prefactor and a sum that is 
essentially a parametrized rank polynomial with only a parameters and one 
variable, sz. We prefer formula ([TJ, though, because we do not want to assume 
invertibility of the b parameters. 

Suppose that the various parameter values a(w) and b(w) are independent 
indetcrminates, and let P denote the ring of polynomials with integer coef- 
ficients in the set of 2 + 4|l / (G)| independent indeterminates {s,z} U {a(w), 
b(w) | to G W(G)}. Let J be the ideal of P generated by the set of 2 \V(G)\ 
products {a(v^)a(v x ), b(v^,)b(v x ) | v G V(G)}, and let n : P — >• P/J be the 
canonical map onto the quotient. Then the only summands of ([1]) that make 
nonzero contributions to itt(M(IA(G)) correspond to subsets T C W(G) with 
the property that \T fl {v$,v x }\ — 1 Vv £ V(G). Each such T is a transver- 
sal of the partition of W(G) into 2-element subsets {v^,,v x } corresponding to 
vertices of G; we denote the collection of all such transversals T(W(G)). Each 
T G T(W(G)) has \T\ = \V(G)\ = r G {W{G)), so s and z have the same expo- 
nent in the corresponding term of itt(M(IA(G)): 



Every generator of J involves the product of two different parameters a(vt), 
b(vi) corresponding to a single vertex v of G. It follows that n is injective when 
restricted to the additive subgroup A of P generated by products 



VteT / \w<£T / 

where k > and T £ T{W(G)). Consequently there is a well-defined isomor- 
phism of abelian groups 7r _1 : n(A) — > A, and we have 



Note that it 1 ttt(M(IA(G)), the image of the parametrized Tutte polyno- 
mial t(M(IA(G)) under the mappings tt and 7r _1 , might also be described as 
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the section of t(M(IA(G)) corresponding to T(W(G)). Either way, formula © 
describes an element of P, where s, z and the various parameter values a(w), 
b(w) are all independent indeterminates. 

Arratia, Bollobas and Sorkin 4 define the two-variable interlace polynomial 
q(G) by the formula 

9 (G) = Yl {x-iyW G)[S]) (y-l)W- rWG)[S]) 

SCV(G) 

f 1 \ \S\-r(A(G)[S]) 

= £ (Hi) <^ 1 ' ISI - 

SCV(G) V 7 

Here r(A(G)[S}) denotes the GF(2)-rank of the principal sub-matrix of A(G) 
involving rows and columns corresponding to vertices from S. 

For T G T(W(G)) let S(T) = {v G 7(G) | u x £ T}; then T h> 5(T) 
defines a bijection from T(W(G)) onto the power-set of V(G). As r G (T) is the 
GF(2)-rank of the matrix 

(columns with u ^ S(T) | columns v x with v G S(T)) 

and the columns are columns of the identity matrix, 

r G (T) = |T/(G)|-|5(T)|+r^(G)[5(T)]). 

It follows that g(G) may be obtained from ir^ 1 itt(M (I A(G)) by setting aiv^) 
I . a(u x ) = a; — 1, &(«</,) = 1, 6(u x ) = 1, s = y — 1 and z = l/(x — 1). These 
assignments are not unique; for instance the values of s and z may be replaced 
by s = (y — l)/<7 and z = cr/(a; — 1) for any invertible a. 



3 Recursive formulas 

In this section we show how a recursive description of the parametrized rank 
polynomial yields the recursive description of the interlace polynomial given by 
Arratia, Bollobas and Sorkin 4 . Recall that a coloop of a matroid is an element 
that is included in every basis, and a loop is an element that is excluded from 
every basis. Suppose r RI is the rank function of a matroid M on the ground 
set W. If w G W then M — w is the matroid on W — {w} whose rank function 
is given by r M (T), for T C If - {w}; M/io is the matroid on W — {w} whose 
rank function is given by r M (T L) {w}) — r M ({w}). 

Parametrized rank polynomials may be calculated recursively as follows: 

1. If is the empty matroid then r(0) = 1. 

2. If M is a matroid on W and w is a coloop of M then r(M) = (a(w) + 
sb(w)) ■ t(M/w). 

3. If w is a loop of M then t(M) = (b(w) + za(w)) ■ t(M — w). 
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4. If w is neither a coloop nor a loop of M then r(M) = b(w)r(M — w) + 
cl(w)t(M/w). 



Suppose G is a graph and v £ V(G) is not isolated; let w be a neighbor of 
v. If v is looped then {a;^ | x G V r (G)} and {v x } U {a;^ | a; ^ v} are both bases 
of M{IA{G)), because I and 



are both of rank |V(G)|, where I' is the identity matrix of order |V(G)| — 1. If 
v is unlooped then {x^ \ x e V(G)} and {v x ,w x } U {x^ \ v ^ x ^ w} are both 
bases of M(IA(Gj), because I and 



are both of rank |y(G)|, where I" is the identity matrix of order |V(G)| — 2. In 
cither case, we see that is not a coloop or a loop of M(IA(G)). 

If v is looped then {v x } U {x^, \ x ^ v } and {w x } U {x^ \ x ^ v} are both 
bases of M(IA(Gj) — v<f,. If v is not looped then {w x } U {x<t> | x ^ v} and 
{v x ,w x } U {x$ | v ^ x ^ w} are both bases of M(IA(G)) — v^. In either case, 
we see that v x is not a coloop or a loop of M(IA(G)) — v$. 

As has only one nonzero entry, the definition of matroid contraction 
mentioned at the beginning of this section tells us that the rank function of 
M(IA(G))/v c f > is the function on W(G) — {v^} defined using the columns of 
the matrix IA(G)' obtained from IA(G) by removing both the column and 
the row corresponding to w. In particular, the rank of the whole matroid is 
|V(G)| — 1. As {xtf, | x ^ v} and {v x } U {x^ \ v ^ x ^ w} are both bases, v x is 
not a coloop or a loop of M(IA(G))/v ( p. 

Let G — v be the graph obtained from G by removing v and all edges incident 
on it. Then IA(G — v) is the matrix obtained by removing the column v x from 



ttt(M(IA(G))) = 

7ra(^)6(« x )r((M(/A(G))/^) - v x ) + ^)a{v x )r{{M{IA{G)) V<j) )/v x ) 
= na(vt)b(v x )r(M(IA(G „))) + nb(v^)a(v x )r((M(IA(G)) V<j> )/v x ). 



7r- 1 7TT(M(7A(G))) = •k- 1 -kt{M{IA{G-v)))+{x-1)-k- 1 -kt{{M{IA{G))-v ! , > )/v x ) 






With the parameter values given at the end of Section 2, this yields 



or equivalcntly, 



q(G) = q(G -v) + (x - 1)^\t{{M{IA{G)) - vt)/v x ). 



(3) 
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Suppose v is looped, and let G" be the graph obtained from G by toggling 
all adjacencies between neighbors of v, and also toggling the loop status of every 
neighbor of v. Consider two matrices 

1 1 0\ /l 0\ 

1 B C and 1 B C\ 
D E) \Q D EJ 

where bold numerals indicate row and column vectors and the overbar indi- 
cates a submatrix in which all entries have been toggled (reversed) . Elementary 
column operations tell us that the two matrices have the same GF(2)-rank. 
Consequently, if v x G T e T(W(G)) then r G {T) = 1 + r G "- v {T - {v x }). The 
definition of matroid contraction tells us that if v x G T £ T(W(G)) then the 
rank of T- {v x } in (M(IA(G)) - v<f,)/v x is r G (T) - 1. It follows that the ranks 
of T - {v x } in (M(IA(G)) - v <jl )/v x and M (IA(G V - v)) are equal. Combining 
this equality with (|3]), we see that if v is looped then 

q{G) = q(G -v) + (x- \)tt- 1 ttt{M{IA{G v - »))) (4) 
= q(G - v) + (x - l)q{G v -v). 

Formula (|4j) is one of the two fundamental recursive formulas for q(G) [4]. 

Deriving the other fundamental recursive formula takes a little more work, 
because the term tt~ 1 itt((M(IA(G)) — v <f>)/ v x) m formula ((3]) does not corre- 
spond to a single interlace polynomial. 

Suppose that v is unlooped and has an unlooped neighbor w. As in [4], say 
that two vertices x, y ^ {v, w} are distinguished by {v, w} if they have different, 
nonempty neighborhoods in {v, w}, and let G vw denote the graph obtained from 
G by toggling all adjacencies between vertices distinguished by {v,w}. 

The matroid (M(IA(G)) - v 4> )/v x has W(G) - v x } = W(G - v) as its 
ground set, and tt~ 1 ttt((M(IA(G)) — v^/v^ includes nonzero contributions 
from elements of T(G — v). Split tt~ 1 itt((M(IA(G)) — v^/v^ into two parts, 
S*0 and S x , with S^, including the contributions from elements of T(G — v) that 
include and S x including the contributions from elements of T(G — v) that 
include w x . 

Consider three matrices 
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Elementary column operations tell us that the first two matrices have the same 
GF(2)-rank, and elementary row operations tell us that the third also has the 
same Gi 7 '(2)-rank. Removing the second row and column from the third matrix 
reduces the GF(2)-rank by 1, clearly. It follows that if v x ,w rj} e T £ T{W(G)) 
then 

r G {T)-l = r G ™- w {T-{wt}). 

Recall that r G (T) - 1 is the rank of T - {v x } in (M(IA(G)) - v^)/v x , and 
the parameter values given in Section 2 include a(v x ) = x — 1 and a(w$) = 1. 
It follows that the contribution of T - {w^} to tt- 1 ttt(M(IA(G vw - w))) is the 
product of x— 1 and the contribution of T— {v x } to tt^ 1 ttt((M(IA(G))—v c /,)/v x ). 
Consequently, if we split q(G vw —w) into two parts, and q x , with q^ including 
the contributions from elements of T(G VW — w) that include and q x including 
the contributions from elements of T(G VW — w) that include v x , then 

(x - i)Sf = q x . 

If ity S T € T(G VW — w), then the corresponding column has only one nonzero 
entry; elementary column operations show that the rank of T in M(IA(G VW — 
w)) is 1 more than the rank of T - {v^} in M (IA(G VW -v — w)). It follows that 

<Z = q(G vw -v-w). 

Formula ((SJ now tells us that 

q(G) = q(G -V) + (X- + S x ) - q(G - v) + q x + (x - 1)S X (5) 

= q{G -v) + q(G vw -w)-q^ + (x- 1)S X 
= q(G -v) + q(G vw -w)- q{G vw - V - w) + (x - l)S x . 



It remains only to discuss S x . Consider two matrices 
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Elementary column operations tell us that the two matrices have the same 
GF(2)-rank. It follows that if v x ,w x e T £ T(W(G)) then 

r G (r)=2 + r G ^ — (T-{v x ,w x }). 

The rank of T - {v x } in (M (IA(G)) - v <t> )/v x is 

r G (T) 1 = 1 + r ™-*-™ (T - {v x , w x }) 
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so considering the parameter values given at the end of Section 2, we see that 
the contribution of T — {v x } to S x is the product of x — 1 and the contribution 
of T - {v x ,w x } to 7t- 1 7tt(M(.L4(G™ ; - v- w)). It follows that 

S x = (x-l)q(G mB -v-w), 

so formula ((5]) tells us that 

q(G) = q{G -v) + q(G vw - to) - g(G™ _«-«;) + (a - l) 2 g(G""" - «- «;). 

This is the second fundamental recursive formula for q(G) given in [4]. 

4 The identity- adjacency- sum matroid 

We call the binary matroid M(IAS(G)) represented by the matrix 

IAS(G) = (I | A(G) | A(G) + /) 

the identity-adjacency-sum matroid of the graph G. The column of A(G) + I 
corresponding to a vertex v is denoted v$, and the other columns are denoted 

00, i> x as in Section 2; then W(G) = {v^, v x , \ v 6 V(G)} is the ground set 
of M(-L4S(G)). Let P' be the ring of polynomials with integer coefficients in the 
set of 2+6 |V(G)| independent indeterminates {s, z}Li{a(w), b(w) \ w € W'(G)}, 
let J' be the ideal of P' generated by the set of 4 |V(G)| products {a(v^)a(v x ), 
a{v<p)a{vTp), a(v x )a(v^), b(v<p)b(v x )b(v^) \ v e V(G)}, and let n' : P' — > P'/J' 
be the canonical map onto the quotient. Then the discussion of Section 2 is 
readily modified to show that (tt')^ 1 tt't(M (IAS(G))) consists only of terms 
associated with 

T'(W'(G)) ={TC W'{G) | \Tn{v lj> ,v x ,v i ,}\ = lVv£V(G). 

The reader familiar with the interlace polynomials introduced by Aigner and 
van der Hoist [T] , Courcelle [T3] , and the author [25] will have no trouble modify- 
ing the discussion of Section 2 to show that appropriate values for s, z and the a 
and b parameters yield all of these interlace polynomials from the parametrized 
rank polynomial r(M (IAS(G))). Notice also that M (IA(G)) is a sub-matroid of 
M(IAS(G), and (tK) - Vt(M (IAS(G))) yields 7r- 1 7 rr(Af(M(G))) by assigning 
a(v^) = and b(v^) = 1; consequently the theory given in Section 2 is contained 
in the one described here. (It is for for expository convenience, not generality, 
that we detail the theory of M(IA{G)) rather than that of M(IAS{G)).) 

5 Comments 

1. As explained by Ellis- Monaghan and Sarmiento [15] . results of Las Vergnas 
p"9l [20] and Martin [21 on circuit partitions of planar 4-regular graphs imply 
that if G happens to be a circle graph obtained from a planar 4-regular graph, 
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then the vertex-nullity interlace polynomial of G may be obtained from the 
"diagonal" Tutte polynomial of an associated checkerboard graph. (The diag- 
onal Tutte polynomial is obtained by setting the two variables of the ordinary 
(non-parametrized) Tutte polynomial equal to each other.) This connection 
cannot extend directly to non-planar graphs, as Martin [2]] p.76] pointed out, 
because the complete graph K 5 has too many Euler circuits to be represented 
in a 5-element matroid. 

Also, Aigner and van der Hoist [T] observed that the vertex- nullity inter- 
lace polynomial of a bipartite graph may be obtained from the diagonal Tutte 
polynomial of an associated binary matroid. This result is connected to the 
preceding paragraph through de Fraysseix's theorem connecting bipartite circle 
graphs to planar graphs [17] • More recently, Brijder and Hoogeboom [11] have 
introduced interlace polynomials for dclta-matroids, and in particular for ma- 
troids. They use them to extend the connection between vertex-nullity interlace 
polynomials and diagonal Tutte polynomials to arbitrary matroids. 

One way to summarize the content of the present note is this: Using pa- 
rameters has the effect of algebraically restricting t(M(IA(G)) to T(W(G)) 
and t(M(IAS(G)) to T'(W'(G)), and unlike restriction to the diagonal Tutte 
polynomial, these restrictions are effective for all graphs and all interlace poly- 
nomials. 

2. The matrix IA(G) appears in pQ and [3J, together with the observation 
that the vertex- nullity interlace polynomial of G is equal to Bouchet's Tutte- 
Martin polynomial of the isotropic system associated with the row space of 
IA{G) 

The matrix IAS(G) also appears in [T], where Aigner and van der Hoist 
showed that their interlace polynomial Q may be obtained by summing over 
submatrices of IAS{G) associated with elements of T'(W'(G)). The content 
of this note came to mind after R. Brijder pointed out the appearance of the 
same submatrices in our work with nonsymmetric modified interlacement ma- 
trices [37] , when we read Bouchet's comment 19! that Eulerian multimatroids 
are "sheltered" by matroids and wondered whether M(IAS(G)) is in general an 
appropriate "sheltering" matroid for the 3-matroid associated with an isotropic 
system with fundamental graph G. Indeed it is! 

3. Some known properties of interlace polynomials can be readily explained 
using known properties of Tutte polynomials. For instance, the analogy be- 
tween pendant-twin reductions for the interlace polynomial and series-parallel 
reductions for the Tutte polynomial noted by Blaser and Hoffman [5], Ellis- 
Monaghan and Sarmiento [15) and the author [251 126) is more than an analogy: 
when v is pendant on w, v x and are parallel in M(IA(G)); and when v and 
w are twins, v x and w x are parallel in M(IA(G))/{v ( j > , u^}. Known properties 
of Tutte polynomials also provide new insights into interlace polynomials; for 
instance, the interlace polynomials of G have activities expansions with respect 
to bases of M(IA{G)) or M{IAS{G)). 

4. In closing: understanding more about the graph-theoretic significance 
of the matroids M{IA(G)) and M(IAS(G)) would help in understanding the 
significance of the interlace polynomials. 
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